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Abstract 

We study two and three jet events with a large rapidity gap at HERA. Unlike in the Ingelman- 
Schlein approach we do not adscribe a structure to the Pomeron. Instead, the coupling of the 
Pomeron to quarks or gluons is taken pointlike, which makes the model easy to test: the only 
degrees of freedom are the coupling constants of the Pomeron to the quarks or the gluons and 
a cutoff procedure to keep the Pomeron-gluon coupling well behaved. 



1 Introduction 



The prediction of the Pomeron as a pseudo particle has inspired a large number of experiments 
in the 60's and 70's. Although effective parameters could be determined under the assumption 
that the particle exists, its existence was never shown unambiguously. Recently the so-called 'large 
rapidity gap' events at HERA Q have caused new excitement about the Pomeron hypothesis. 
These events show jet activity far away from the beam pipe, but there is no trace of a color string 
between these jets and the remnants of the proton. Actually it is not excluded that in many or 
even most cases the proton remains intact. Hence the events carry the signature of a diffr active 
mechanism in which a color neutral object reacts with the virtual photon that comes from the 
electron. The outcoming events can be classified by the number of jets they show. This is similar 
to e~e + — > hadrons with which these events can be shown to have much in common. We will 
comment on this in the results section. 

Thus far most theoretical analysis has concerned the two jet events. Various models have been 
tried, ranging from a full calculation of a gluon pair exchange @ and [|, to a phenomenological 
description of Pomeron structure functions [|| ||] and even to a gluon photon fusion scheme in 
which the resulting quark pair somehow sheds its color and leaves it behind inside the proton 
fragments l|. In this paper we will produce the Pomeron with a classically parameterized diffractive 
mechanism [pj. Having produced it, we will couple it to quarks and/or gluons with a pointlike 
coupling. This has been tried before for two jet events, but then in the case of a small Q 2 of the 
photon ||; see also ref ||. In this paper we are however interested in the region of large Q 2 in 
which the outgoing electron gets observed clearly inside the detector. In addition we will have a 
look at three jet events. The aim of the paper is then to see how well the events can be used to 
determine the coupling of the Pomeron to the quarks (and/or gluons). 

There are various motivations for this model. First it is the simplest model one can think of, 
and hence it is necessary to have at least a look at it. Another motivation has to do with the fact 
that the model contains explicit calculations for three-jet production allowing us to make simple 
predictions for broadening effects in the multihadronic diffractive final states. In addition the 
matrix elements can be implemented in a Monte Carlo programme using the Lund fragmentation 
scheme. Thus, our model can be directly compared to experimental data. In fact, as it will be 
shown in a forthcoming publication fi~0|| , the model gives a good description of the data. 

2 The computation 

Considering the fact that we insist on seeing the electron, we can treat it in the normal way. This 
means that it is just another particle in the trace calculations, and we do not have to take extra 
precautions with respect to numerical stability of either the matrix element or the kinematics. 
Because of this there is absolutely no need to use an equivalent photon approximation. If the 
electron were to stay in the beampipe the situation would be different. In the forward direction 
terms of order m 2 /s are important and seemingly leading terms cancel each other to a high degree 
of accuracy. This not being the case we can work out the trace of the Feynman diagrams by brute 
force using the program FORM |]llj and evaluate the resulting FORTRAN code without much 
worry. The two jet diagrams of figure 1 are easy to compute. 
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Figure 1: The 2 diagrams for Pomeron induced two jet production. The 
Pomeron is denoted by the zigzag line at the left bottom of each diagram 

In this case there is only one 'new' parameter. The strength of the coupling of the Pomeron to 
the quarks: gp qq . Basically the normalization of the reaction fixes this constant. The cross section 
is easily calculated. If one were to ignore the polarization of the photon and the 'mass' of the 
Pomeron the matrix element of the subprocess 

7(<Z) + P(p) - + q(k 2 ) 

(in brackets the momenta of the particles) then its cross section, denoted by <t 7 p would be 

2 9 2 Pqq a s-q 2 x , ^d 3 h d 3 k 2 

da lP = —5{p + q-h- k 2 )—— — — . 

7T Ul Zfcio 2k 2 

Here t,u are the standard Mandelstam variables (t is defined by P and k±), and g Pqq is defined by 
the effective interaction Lagrangian C c r = g Pqq qqP- In practice we did not ignore the polarization 
of the photon and the 'mass' of the Pomeron. Hence the formula for the full process is a little 
longer: 

do P ~ d ^ d ^ d "Pe Ut 5Hv p , , out )( 3 2 (Mj- g 2 ) 2 + P 2 2 2 2 1 1 

d<JeP ~ 2E 1 2E 2 2E°"t d [Pe + F kl k2 Pe )[ q ut + {Q M n 2 + t 2j 

,00 „ A , Iw Pe-fcl (P 2 ~U) Pe -k 2 (P 2 ~t) 

+ ^' M u + t )[ t + u } 

M 2 11 q 2 Q 2 

-8(P . „.)._£ _ 8P 2 pe . hi pe . h _ ( _ + _ ?) _^__ (1) 

We use here Mx for the diquark invariant mass, p e for the incoming electron and p° ut for the 
outgoing electron. 

The three parton diagrams, responsible mostly for the three jet events, are given in figure 2. 
Here the formulae are much more complicated, but with the use of FORM we obtain an expression 
of about 800 terms which can be written into a file in such a way that it can be included into a 
FORTRAN program directly. We have now another 'new' parameter: g Pgg , the coupling of the 
Pomeron to two gluons. Its size can be fixed by looking at the ratio between two and three jet 
events. It should be noted that because the Pomeron is taken to be a scalar particle and because 
we assumed that the quarks are massless, the six quark-coupling induced diagrams do not interfere 
with the two gluon-coupling induced diagrams. 
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Figure 2: The 6+2 diagrams for Pomeron induced three jet production 

When the Pomeron is coupled to two gluons we use a gp gg F llv F^y coupling. This presents 
some problems. Because this coupling causes an unphysical high energy behaviour, due to the 
extra dimension, this is 'repaired' by a form factor. For this form factor we took the expression 
1/(M 2 — t^) 1 ^ 2 in which t% is the invariant mass of the internal gluon. It has a negative value. M 2 
is a free parameter which is expected to be of the order of lGeV 2 . Hence we used this value. The 
presence of the formfactor is noticeable in the differential cross section with respect to Mx which 
is the invariant mass of the three jet system. In the case that we do not use the form factor we 
get a clearly unphysical behaviour for large values of Mx- It is roughly constant until phase space 
cuts it off. With the formfactor a more physical behaviour is obtained and the mass distribution 
looks similar to the corresponding one for the 6 diagrams with the quark coupling. 

For the production of the Pomeron we used a phenomenological fit that has been obtained by 
Donnachie and Landshoff Q. This way the Pomeron flux in a proton can be written as: 

fp(x P ,t) = -^{Fxit)) xp (2) 

with 

4m 2 , -At to 

A = 2.8 

B = 0.7 

Pa = 3.24GeV~ 2 

a(t) = l + e + a't 

e = 0.085 

a = 0.25 (3) 

It is of course rather easy to implement such a flux factor, provided that the integrations over phase 
space are done numerically. It is equally possible to change the intercept e to a value close to 0.2 
as suggested by recent HERA data. 
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Next we need to do the integration. This is done by means of a Monte Carlo program VE- 
GAS [O] . It is a so-called self learning program that integrates over a unit cube of a given number 
of dimensions. Hence we have to transform points in such a cube into points in phase space. The 
routines that do this are called the kinematics routines. In the current calculation we use some 
'standard' two and three body decay routines. Both routines assume that the 'decay' takes place 
in the CM frame of the decaying object. This means that at various points we need some trivial 
Lorenz transformations. The two body decay routine is used for the reaction e~P — > e~ X in which 
P is the Pomeron and X represents the 2 or 3 jet system. It is also used for the subreaction 
7*P — ► 2 jets. The three body decay routine first selects two Dalitz variables and then supplements 
them with the proper CM angles. It is used for the subreaction 7*P — * 3 jets. Because everything 
happens at large angles such routines are sufficient. 

Of course integration would be meaningless if we would not apply acceptance cuts. One has to 
worry about two sets of cuts. In the first set we try, at the parton level, to imitate the experimental 
analysis. In this case we do as if the partons are the jets of the experiment and we copy the analysis 
cuts of the experiment. In the second case we prepare the program for event generation allowing the 
experimental analysis programs to apply fragmentation into hadrons. We can then do the analysis 
at this level. Alternatively the resulting particles can be sent to the detector emulator programs. 
This last step would be part of a full experimental analysis. In the case of hadronization the cuts 
have to be stricter, because the fragmentation will make things less clear. There will be events that 
pass the analysis cuts this way that did not pass them in the theoretical analysis, et vice versa. In 
general the last method is superior, but it requires the application of the full experimental detector 
simulators. It also requires many more computer resources. The first method however is usually 
sufficient to give a good idea of what is happening, how distributions will differ between various 
models, etc. The study of the hadronization effects, without using the detector emulators, gives 
already a good indication of whether things will change much in a full analysis. The cuts we use 
in the calculations at the parton level are: 

• Rapitity(r?) of the scattered electron: > —3.7. 

• Energy of the scattered electron: > 5GeV. 

• Q 2 of the electron: > 3GeV 2 . 

• Mass of the outcoming hadron system(Mx): > 3GeV. 

• Angle of the hadron system: cos 9 < 0.9. 

• Rapidity of the individual partons: < 3. 

• Transverse energy of the individual partons: > 0.2GeV. 
. M]jlM\ > Y min = 0.05. 

• Beam pipe cut: Each parton has 9 > O.lrad. 

For the calculation that involves the hadronization we used 

• For the jets apply the cuts for the individual partons. 

• Only particles with 9 > O.lrad are used. 

• There should be at least 4 particles. 
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• All particles should have a rapidity < 2.5. 

• Angle of the hadron system: cos 6 < 0.75. 

• Q 2 of the electron: > lOGeV 2 . 

• Mass of the outcoming hadron system(Mx): > 5GeV. 

The hadronization was performed with the program JETSET 7.4 [|l3|]. The parton flavors were 
selected to be d:u:s = 1:4:1 and the decays of some of the hadrons (mainly the strange ones) were 
suppressed. 

Because we use a Monte Carlo integration method, it is easy to make a large variety of his- 
tograms. In addition we can use an additional set of programs after the integration that makes 
the program behave like an event generator in which the events all have weight one. Hence the 
programs are suitable for a detailed experimental analysis. The number for the total cross section 
can usually be obtained with an accuracy of a few tenth of a percent during the integration phase. 

3 Results 

In diffractive e — p events, the absence of the proton remnant considerably simplifies the analysis 
of the multihadronic final state. If interpreted in terms of 7* — P collisions, the Pomeron fractional 
momentum £, with respect to the proton beam, can be easily obtained from the relation 

? M\ + W 2 1 ' 

where Mf-, Q 2 and W 2 stand for the usual DIS kinematic variables. 

We propose to establish an analogy between 7* — P scattering at a centre of mass given by Mx 
and e + e~ annihilation at a centre of mass energy yfs. In this section, we study the dependence of 
shape variables describing the multihadronic final states as a function of Mx, in order to asses the 
importance of the two different broadening mechanisms in the model 

ft gluon radiation 

4k Pomeron splitting into two gluons 

To 0(a s ), the cross-sections for the three types of processes considered in this paper, within 
the cuts described in the previous section and for a fixed value of the strong coupling constant 
a s = 0.25, are: 



quark pairs 



(43.83 ± 0.08) • g p nb for two-jet production due to the direct coupling of the Pomeron to 



• (3.240 ± 0.003) • 9p nb f° r three-jet production from gluon bremstrahlung to the procees 
above 

• (0.2705 ± 0.0009) • 9 Pgg nb for three-jet production from the Pomeron coupling to gluon pairs 

In Fig. 3a we show ~ d d ^ x for these three reactions. Notice that three-jet processes through 
the Pomeron coupling to gluons tend to give rise to larger masses then those coming from brem- 
strahlung. In Fig. 3b we show the thrust [14| distribution for three parton final states in the 7 — P 



centre of mass frame. Again, the tail in the thrust distribution associated with final states due to 
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Pomeron coupling to two gluons, dashed line, is harder than that associated to bremstrahlung from 
the Pomeron to two quarks coupling, solid line. In fact, the latter is in rather good agreement with 
the analytical expression derived by Ellis, Gaillard and Ross for e + e~ annihilation 15]. 




Xq Xg 

Figure 3: The distributions at the parton level. 



The solid line refers to the production of three partons, using the Pqq vertex. The dashed 
line corresponds to the Pgg vertex. The dotted line in the figure at the top left refers to 
the production of two partons. 

To deepen our insight into the roles played by gluons in the two different three-parton processes 
we are considering, it is instructive to look at the parton fractional momenta, defined as usual by 
x q = 2E q /T,Ehadr and x g = 2E g /T,Ehadr- They are presented in Figures 3c and 3d respectively. It 
is striking to see that while bremstrahlung gluons are soft partons as expected, those coming from 
Pomeron splitting are much harder. 

In order to study to what extent the properties exhibited above are masked by fragmentation 
effects, we have calculated the thrust distributions in the hadronic centre of mass, after fragmen- 
tation, for various ranges in Mj. 
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Figure 4: The thrust distributions at the hadron level. 

The solid line refers to the production of three partons, using the Pqq vertex. The dashed 
line corresponds to the Pgg vertex. 

The results are shown in Fig. 4. Notice that while at low masses, both thrust distributions, i.e. 
the one stemming from the P — > qq coupling and that due to P — > g g are similar, at high masses, 
the latter gives rise to multihadronic final states which are more spherical than the former. 
To summarize and conclude. We have calculated the cross-sections to 0(a s ) for diffractive e — p 
scattering mediated by Pomeron emission. We have considered a direct coupling of the Pomeron 
to quark and gluon pairs. If the Pomeron couples preferentially to gluons as recently indicated 
by experimental data |l6|| , the multihadronic final states should exhibit, as the hadronic masses 
become large, broadening effects much stronger than those associated to gluon bremstrahlung in 
e + e~ annihilation [17]. 



One of the authors (J.V.) likes to thank the Universidad Autonoma in Madrid for its kind 
hospitality. 
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